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Abstract 

A new semiclassical approach to ionization by an oscillating field is presented. 
For a (5-function atom, an asymptotic analysis is performed with respect to a 
quantity h, defined as the ratio of photon energy to ponderomotive energy. This 
h appears formally equivalent to Planck's constant in a suitably transformed 
Schrodinger equation and allows semiclassical methods to be applicable. Sys- 
tematically, a picture of tunneling wave packets in complex time is developped, 
which by interference account for the typical ponderomotive features of ioniza- 
tion curves. These analytical results are then compared to numerical simulations 
[H and are shown to be in good agreement. 
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1 Introduction 

In recent times, there has been a lot of effort dedicated to a better understanding of 
ionization by strong laser fields (reviews e.g. in [@, §]), especially since the discovery 
of nonperturbative effects like above-threshold ionization ATI (for a review see [Q), 
the sensitivity of ionization rates, and stabilization in superintense fields 0, ^, 0- But 
there succeeded no analytical solution, not even to the simplest model, i.e. an electron 
bound by an attractive 5-potential in the presence of an oscillating electric field. One 
of the main reasons is that there exist two separate regions. Inside the atomic core, 
the binding potential is predominant, and outside, it is vice versa the electric field 
that dominates (this is also the main point that makes perturbation theory work so 
poorly). For both regions, the distinct propagators are known exactly, but they cannot 
be combined to solve the ionization problem exactly. 

An appropriate way of addressing this problem nonperturbatively is the semiclassical 
method, which we will use to construct the semiclassical propagator for ionization of 
such a 5-function atom. The choice of this model has three advantages: first, there 
is a clear distinction between inside and outside the atom, so that there exists no 
intermediate region. Second, there is only one bound state for the ^-potential, so there 
arise no difficulties with intermediate resonances and induced resonances by AC-Stark 
shift as it happens in real atoms. In this point, our model resembles photodetachment 
of a negative ion especially H". Third, the exact problem can be reduced to a 
Volterra type integral equation in time, for which accurate numerical solutions can be 
computed [0] and allow precise tests. 

The above described problem (as well as its more general settings in three dimensions 
with more realistic binding potentials) has been treated in the literature in several ways. 
The so-called Keldysh-Faisal-Reiss (KFR) approach |P, [l^, consists in expressing 



the exact propagator If in terms of the known Volkov propagator f/^ [|I^ (for a free 
electron in the electromagnetic field) and of Vs, the atomic binding potential: U = f/^ + 
f/^ Vs U . The unknown U on the right hand side is approximated by f/^ (equivalent to 
Born approximation), and matrix-elements for ionization are calculated between the 
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ground state and so-called Volkov-states [|l^ in the continuum. This approach has 
been extended and refined by many authors, and we will compare our results with two 
such typical extensions [|T^ |T^. 

Another approach, the so-called two-step model [|15|, |16], |l7l, clearly distinguishes be- 
tween ionization first and classical propagation in the laser field afterwards. This 
proved to be very useful especially in calculating high-harmonic generation |T^. This 
separation into two steps will be used in the following, but now justified in a fully 
semiclassical context. 

In addition, there exist several other approaches; a very common method is using 
Floquet theory [l^, which explicitly incorporates the periodicity of the time- 

dependent Hamiltonian. 

Our issue here is not to obtain better results for a simple model, but to gain better 
physical insight into the mechanisms of ionization processes using semiclassical meth- 
ods. 

This paper is organized as follows: After basic definitions, we demonstrate characteris- 
tic elements of typical ionization curves which we want to understand semiclassically. 
Using a "sum over classical paths" technique, the total semiclassical propagator is con- 
structed successively by identifying the paths which are relevant for ionization. First 
for the part remaining bound, and then for the free wave packets stemming from time- 
dependent tunneling in complex time. From this semiclassical propagator, the total 
ionization rate is derived and compared to numerical simulations as well as to other 
theories in the literature. 
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2 Basics 

2.1 Definition of the model 

We want to study the ionization of a one-dimensional 5-function atom Vs = —<yS{x), 
which possesses exactly one bound state with energy Eq = — (for the three- 



Vo = —fix cos{ujt) in the so-called dipole-approximation. The parameter a stands for 
the strength of the binding potential, /i for the amplitude of the applied electric field, 
and uj for its angular frequency. The Schrodinger equation in atomic units {h = m = 1) 
is 



From these three parameters, there can be derived (due to the scaling properties of 
(|l|)) two independent, meaningful, and dimensionless quantities: 

first, there is z, the ratio of ponderomotive energy Upond = A*^/(4t^^) to photon 
energy u: z = fi"^/ (Auj^). Upond is the mean kinetic energy of a free electron in an 
oscillating field. 

and second the so-called Keldysh-factor [0 7 = au/^, which is the ratio of the 
(adiabatic) tunneling time to the period of the applied field. This 7 characterizes 
the ionization process; 7 -C 1 corresponds to (adiabatic) tunnehng, and 7 ^ 1 is 
better described in a pure multiphoton frame |^ p^ . 

These are the two quantities by which the model will be described below. A third, 
but no more independent quantity is njo, the ratio of necessary ionization energy Eq 
to photon energy: 



dimensional analog cf. |22]). This atom is exposed to an oscillating electric field 




(1) 



(2) 
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2.2 Transformation of the Schrodinger equation 

Using the following coordinate and time transformation, we cast the above Schrodinger 
equation in a very suitable form: 

x' = — X , t' = ut (3) 
/i 

Using the scaling relation 6{ax) = 6{x)/\a\, which is identical to the scaling behaviour 

of the Coulomb-potential, and omitting the primes, we arrive at: 

d / 1 d'^ \ 
ih—ip=i--h^—-h'j5{x)-xcos{t)j'4j (4) 

Here h := uj^/fi"^ = 1/(4 z) is written very suggestively to indicate that we have obtained 
a parameter h formally identical to Planck's constant h in ordinary quantum mechanics 
using SI units. Of course we can give h any value we like. Restricting ourselves to strong 
fields where Upond ^ ^, i-e. z ^ 1, h can get arbitrarily small. This will allow us to use 
the normal semiclassical methods, exploiting the formal analogy between the parameter 
h, introduced above, and Planck's constant h. Still remarkable is the appearance of 
the factor h in front of the binding 5-potential, its implications in the semiclassical 
limit will be derived below. 

The (normalized) ground state wavefunction for the (5-function atom without applied 
external field is 

Mx) = -y/^exp(-^|x|) (5) 
Hoipo = -yV'o (6) 
i'oix,t) = ipo{x) exp{i—t) , 

which leads to the stronger a localization the smaller h is. In the limiting case h ^ 0, 
ipo^x) approaches (appropriately scaled) the spatial Dirac ^-function: 

I^HMx) - S{x) (7) 

In the following, this approximation will be used only for calculation of scalar products, 
so that no mathematical ambiguities will arise. 



3 Prominent Semiclassical Features in Ionization 
Rates 
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Figure 1: Numerical ionization rates, raw (dotted) and smoothed (thicker) versus 
number of ponderomotive photons compared with WKB-background (thick, smooth 
decaying curve). The Keldysh parameter 7 and not the depth of the binding potential 
is kept fixed. 



The ionization rates from a numerical analysis [|T| of the model show certain charac- 
teristic features. In figure |T], we see the raw results (dotted), these results smoothed 



(thicker) and the adiabatically averaged WKB-value (see eq. ([T6|) , depicted thick, 
monotonically decreasing). The ionization rate F is shown against the ponderomo- 
tive energy (over photon energy). The WKB-rate accounts well for the background, 
but on the actual rate, there is a periodic modulation and a lot of fine-structure super- 
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posed. The (slow) modulation becomes obvious after smoothing the raw data (using 
a Savitzky-Golay filter technique [^) in order to eliminate the fine structure and to 
work out this fundamental modulation. 

The period of this modulation can be understood by so-called channel closing ar- 



guments |26|, but there is no argument for the amplitude of modulation. The fc-th 
channel means a multiphoton ionization by k photons with energy balance (here in 
atomic units) 

ku = UioUJ + ZUJ + Ekin (8) 

riio ijj is the energy required for ionization, z w is the (ponderomotive) energy U^ond con- 
tained in the electron oscillation in the electric field, and is the additional kinetic 
energy the free electron gains in this ionization process. This channel is energetically 
only allowed for > 0, otherwise it is forbidden. 

The threshold for the fc-th channel is defined using the condition E^^in = 0. In terms 
of 7 and z, equation (|) at threshold is 

k = + Zk (9) 

and the specific value z^ at threshold becomes 



k_ 

1 + 272 



yielding a period Az = 1/(1 + 27^) in z, i.e. the period of approximately one half in 
the example of figure |l|. 

An important point to note is that the numerical results give strong evidence for regular 
behaviour at threshold, whereas the usual prediction of appropriate theories (e.g. those 
of KFR type) is a divergence at threshold. 

Here it is important to note that the features described above are definitely not re- 
stricted only to the 5-function atom. Numerical simulations for various model po- 
tentials in the literature exhibit similar features, often with remarkable quantitative 
correspondence (see figure 8 in |1[], comparing the (5-function atom to a smoothed 
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binding potential V{x) = — exp{—\x\) / y x'^ + Xq as used by Greenwood and Eberly 

im)- 

In the following, a semiclassical theory is derived that accounts for the information 
contained in smoothed rates. The WKB-background as well as the periodicity, ampli- 
tude, and phase of modulation are contained in a single, divergence free theory, which 
constructs the propagator using the semiclassical sum over paths, cf. e.g. . 



4 Quasi— Energies 

4.1 Calculation of WKB-coefRcient 

The main effect of applying an external field to an atom is that the bound state becomes 
metastable and tunneling can occur. In the static case, this tunneling rate D can be 
approximated using the normal WKB-coefficient for the corresponding barrier. For a 
linear potential barrier V{x) = —rjx and a given (negative) energy Eq, the well-known 
expression for D is 

-Eo/ri 

h 





£)(77) ^ exp ( — / \p{x)\dx] (11) 



Using the (imaginary) local momentum p{x) = \J2{Eq + rjx) and the ground state 
energy Eq = —7^/2, this evaluates to 

Dir])^exp(-l^) (12) 



3T]h 

Of course, this approach gives only the exponential part, but since the preexponential 
part is well-known from the literature for this simple case (e.g. [|14|, |2^), we can adopt 



it from there. This factor is just twice the atomic frequency, in our units 2Eo/h 

m = t^^p(-H) (13) 



h \ "irjh 

If we consider a time-dependent external field, the parameter 77 becomes time-dependent 
too and represents the instantaneous strength of the applied field at the moment the 
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ionization takes place: rj = \ cos(t)|. If the tunneling process occurs on a much shorter 
time scale than the period 27r of the oscillation, it is a good idea to consider the 
ionization taking place adiabatically, i.e. calculate the instantaneous ionization rate 
D{\ cos(t)\) and average it over a whole period. This case corresponds to 7 <^ 1, i.e. 
the Keldysh factor characterizing the ionization is quite small. 

In calculating the cycle-average D over a cosine period, one has to integrate and to 
normalize subsequently 

D = — /^expf--- — -CTTT)dt (14) 
27ri /i 3\cos{t)\hJ ^ ' 

Because we want to examine the asymptotic case /i — > 0, we best evaluate this integral 
using the method of steepest descent (also called saddle-point integration, cf. appendix 
0) with respect to h . The derivation of the exponential argument with respect to time 
t gives 

d ( 2 73 \ _ 23 sin(t) 



dt\ 3cos(t); 3 cos(t)^ 
The relevant times are the zeros of this expression, i.e. all multiples of tt. These are the 
instants where the electric field strength | cos(t)| is a maximum. For symmetry reasons, 
all these instants are equivalent and it is sufficient to evaluate the above integral around 
one such instant using the method of steepest descent. The result is 



This shows that the average D is the instantaneous ionization rate at maximum field 
strength, up to a preexponential factor. 

Note that the method of steepest descent becomes exact in the limit of vanishing h. 
This allows a very interesting interpretation. In this case, the ionization effectively 
takes place only in the vicinity of the instants of maximal field strength rj = 1. This 
means that there exist ionization bursts, separated by half a period, between which 
practically no other ionization occurs. In the following, this property will be used to 
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construct a scenario of propagating wave packets, emerging from these peak instants, 
freely propagating afterwards, and interfering with one another and with the part of 
the wave function remaining bound. 



4.2 AC-Stark effect 

The second effect of applying an external field is the so-called Stark effect, in the 
oscillating time-dependent case called AC-Stark effect . We will treat this effect adia- 
batically as well (cf. ||30|), the well-known value (e.g. |jl^ ^ ^) for the instantaneous 



energy-shift E^'^' is E^^^ijii) = — 5/1^77^/(87^). Cycle-averaging results in 

2 ^' ' I674 ^ ' 

This effect means an additional phase factor in the propagator exp{—iHt/h), whereas 
the tunneling rate, calculated by taking absolute squares of the wave function, is not 
affected by E'^'^. Note that the influence of (|1^) will be quite small in the following 
because of its proportionality to h"^. 

In order to express tunneling and Stark shift together, it is useful to write the exponen- 
tial decay of the bound state using an imaginary contribution to the total energy 
E™. Setting 

E\v) = -^t^w(-l^)=-4Div) (19) 



2 "V 3hr]J 2 
E^{r]) = Eo + E^''{r])+E'{r]) , (20) 

the adiabatic development of the ground state can be described by the propagator 

U\tj) = exp(--^|E"^(|cos(t)|)^/t) Po, (21) 



using the total quasi-energy E"^ and the projection operator Pq, projecting onto the 
ground state. This adiabatic description is useful if we want to describe the propagation 
of the wavefunction for arbitrary times. If we restrict ourselves to considering only full 
cycles, we can use the appropriate averages. 
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= Eq + E^^ + E^ (23) 
u\tj) = exp(^-^E"^ty)Po (24) 

Applying this propagator to the ground state and taking absolute squares results just 
in the exponential decay with rate D. This will be sufficient for the forthcoming 
considerations. For tj = 2 kn, k = 0,1,2, both propagators are of course identical, 
due to the very construction of the average. 

When comparing these expressions with the numerical results described later, we will 
see that they can account for the monotonic background of the ionization rate (see 
figure [^), but if we want to explain the quasi-periodic modulations, responsible for the 
nonmonotonicity of the rate, we have to go further in our semiclassical description. 



5 Semiclassical Propagators 
5.1 General construction 

To construct the semiclassical propagator outside the binding potential, we start with 
the (formal) path integral expression 

*/ 

U{x,tf-y,t,) = |Px(t) exp Qs[x(t)]) (25) 

u 
h 

S[x{t)] = J L(^x{t),x{t)^dt (26) 

u 

L = T-V = Lo-Vs = T-Vo-Vs (27) 

Here T is the kinetic energy operator, L is the full Lagrangian, and Lq is the Lagrangian 
for the electric field Vq alone, with the binding potential Vs excluded. The usual pro- 
cedure in the semiclassical limit /i is to find the stationary paths with 6S = 0, 
yielding the classical paths by means of the Euler-Lagrange equation |3^. The remark- 
able point here is that Vs = —h'y6 contains a factor h. This becomes important for the 
semiclassical limit, because this h cancels in the exponent iS/h, and consequently this 
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part of the phase does no more fluctuate arbitrarily fast for non-stationary paths in 
the semi classical limit. 

Applying the analogy to saddle-point integration in function space (cf. appendix 
we notice that we only have to vary Sq = J Lq in order to find the stationary paths to 
6Sq = 0. This condition gives the classical path Xci{t) to Lq by means of the Euler- 
Lagrange equation 

subject to the boundary conditions imposed by the path integral. 

XciiU) = y , Xci{tf) = X (29) 

Xd{t) = Xci{t\x,tf;y,ti) (30) 

X — y + cos(t f) — cos(ti) , , , , 

= - cos(t) + cos(ti) + 1/ + ^ -^if ^(t-ti) (31) 

tf - ti 

Vs in the full Lagrangian L only accounts for an additional phase factor exp(i70) to 
the propagator 



6ixciit))dt (32) 

u 

= EV|icK^^o)h (33) 

'■0 

where the t^ denote the zeros of the classical path Xditl) = 0. This phase (p jumps 
every time the classical path Xciit\x,tf;y,ti) crosses the 5-potential at the origin. 

The result for the semiclassical propagator U^'^ is 



1 92 i 

"'°^'''^y''^^ = Vl^f exp(-S„)exp(n0) (34) 

1 ■ 

exp(- / LQ{xci{t),Xd{t))dt]ex-^{i'-i(t)) (35) 



l2mh{tf - 1,) V/i 
The same result is derived in appendix |B| using the WKB-ansatz. 
In general, one would have to include so-called Maslov phase factors |Q, but we can 



omit them because we do not encounter any caustics in this problem. Since Vq 
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— xcos(t) is linear in x, the appropriate semiclassical propagator for Lq is identical 
to the exact one, namely the well-known Volkov propagator ||12[. The result can 
now be understood as the Volkov propagator plus additional phase jumps for every 
crossing of the origin. 

W = exp(i70) (36) 



5.2 Special tunneling propagator 

There is one important point to note; the description of the propagator using regular 
classical paths is only justified after the electron has tunneled out. So in order to 
describe the tunneling paths, which classically do not exist, one has to modify the 
above description; a common method is to introduce complex time and coordinates 
(cf. appendix ^ and P5| , 



In our case, we know the bound state i/jq = y7//iexp(— 7|x|//i), which formally equals 
a free wave exp{ipx/h) with complex momentum po = ±27. This is consistent with a 
negative energy E = Pq/2 = —7^/2, which is just the ground state energy Eq of the 
(5-potential. 

We choose the following complex boundary conditions (a similar reasoning appeared 
in 1^) for the complex tunneling path Xt'- 

the initial momentum (imaginary part Im considered only) 

Im(xT(to)) =Po = +17 (37) 

(plus sign chosen to ensure exponential decay of wavefunctions and not growth), 
and the initial position 

xrito) = y (38) 
Additionally, xt must fulfill the final condition 

xritf) = X , (39) 
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which is the boundary condition at the end of the path. The additional free constant to 
is necessary because we impose three boundary conditions. But the ordinary differential 
equation (^), which xt must obey, is of order two, and therefore has only two free 
constants. 

The idea is that tunneling takes place in the imaginary part between t = to and t = tj, 
and free propagation U^'^ (under the influence of the oscillating electric field) takes place 
vice versa in the real part. This interpretation is allowed by the usual decomposition 
rules for semiclassical operators (cf. |^8|, ^). 

The general expression for such a path xt fulfilling the (now complex-valued) equation 



of motion (pSf) is 



xrit) = - cos{t) + cos{to) + y + vo{t - to) (40) 

xrit) = sm{t)+vo (41) 

The first condition lm(xT(to)) = Po yields 

i7 = Im(^sin(to) + wo) (42) 

Now we see the meaning of to] it must account for the complex boundary condition 
and so we set 

to = i arcsinh(7) , (43) 



in order to allow vo to remain real (see also fS^] for this result). The second boundary 
condition (|38| ) is fulfilled trivially by the ansatz (0), and from the third conditon we 
obtain 

_x-y + cos{tf)-cos{to) 

Vq - 7 (,44J 

tf - to 

Using the equations (^,|^,0), we can easily construct the complete propagator U'^ 
containing tunneling as well as propagation in the electric field. The result is just the 
analytic continuation of our former result U^'^ in eq. ( ^5]) 

1 i *^ 

U'^{x,tf; y, to) = ^ = exp (^- / Lo{xT{t),XT{t)) dt^ exp(i70[xT]) (45) 



^2mh{tf-to) 



to 
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Note that the argument of the square root in the denominator is now truely complex, so 
that we have an ambiguity in choosing a certain sheat of the complex root. We decide 
to define ^ r exp{i{p) = —^/rexp{i{p/2), ip G [0,27r]. Note further that the tunneling 
starts at to = arcsinh(7), and in order to reach to, we first have to propagate the 
ground state from t = to to using U (to), the analytic continuation in time of (p^. 
This path in the complex t-plane is depicted in figure 0. 

According to the composition rule for propagators, the complete propagator If^ for the 
ground state from t = to t/ is 

U%tf,0) = U^{tf,to)lf\to) (46) 
or in coordinate representation 

+ 00 

If'ix^tf, y,0)= J U^{x,tf- z, to) u\z, to; y, 0) dz (47) 

— oo 

5.3 Generalization to other ionization bursts 

The result of the previous section can be easily generalized to later ionization bursts. 
Figure |^ shows in the complex t-plane that the electron does not become instanta- 
neously free, but rather propagates from t = to t = to according to the propagator 
U (to) (valid inside the binding potential), and then tunnels and propagates from to 
to tj = 27r according to the complex-valued propagator If'^ valid outside. The above 
line integrals over the depicted paths are path-independent so there does not ex- 
ist a unique path. The important point is the propagation with different propagators 
from distinct starting points to distinct endpoints. This describes the first ionization 
burst, but in order to describe the wave packets emerging from the bursts at times 
tk = k7i,k = 0,1, 2, 3, kn <tf, one can repeat the above calculations. 

The propagator U^itf, 0) for the wave packet stemming from t = fcvr is just 

U^{tf, 0) = f/^(t/, kn + to) u\kn + to) (48) 

In this notation, Uq is identical to the above f/^. The interpretation is that the electron 
remains trapped from t = tot = kn and then propagates into the complex t-plane up 
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Figure 2: Path in tlie complex t-plane, describing the evolution of the first wave packet 
stemming from t = 0. Since the line integral for the calculation of the semiclassical 
propagator is path-independent, these paths can be chosen at random. Important 
is that the propagation from t = to t = to uses another propagator U than the 
propagator f/^ afterwards to tf = 2 n. 

to t = kn + tQ, both according to f/ . It then tunnels and propagates from t = kn + Iq 
up to t = t/, according to (^51), i.e. the (complexified) Volkov propagator plus 
phase jumps. Figure ^ contains the four paths corresponding to four wave packets 
created at t = 0, vr, 2 vr, 3 vr, which interfere at = 4 vr. Once again, these line integrals 
are path-independent so that one can choose these paths at random. 
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Figure 3: The various patlis in the complex t-plane for four wave-packets are depicted. 
The path for the wave packet IT goes from to k tc + to and afterwards 

to tf = Att. Again we have path-independence for the two parts. 

6 Interference between Paths 

The electron has two possibilities, it can tunnel and propagate as well as remain bound 
by the binding potential. The quantum mechanical amplitudes for both processes are 
known and they can be added in order to obtain a better description of the evolution 
of the system. This is just the "semiclassical sum over classical paths" method, e.g. 

The full propagator U is the sum of the propagator U , valid for the part bound by 
the (5— potential, and the 's, the propagators for tunneling at t = kn and (free) 
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propagation afterwards. 



U{tf,0) = u\tf)+ U^,{tf,0) (49) 

A:=0,l,2,.-- 



6.1 First period 

For simplicity and notational reasons, we will first consider only the wave packet origi- 
nating from the ionization burst at t = 0. The other propagators can be added as well 
and the resulting integrals can be evaluated using the same techniques as described 
below. We will restrict ourselves to examining the wave function after full periods 
tf = 2n7r, and we will deal mostly with just one period. The second burst occuring 
at t = vr during the first period is of secondary importance, because this free electron 
follows a classical trajectory Xd = a + b{t — n) — 1 — cos{t) and its center is about —2 
to the left at tf = 2mx. Therefore the overlap of this wave packet with the ground 
state ipQ can be neglected. 

This effect as well as the influence of considering several periods and wave packets 
will be demonstrated, when evaluating the analytic expressions derived below and 
comparing them to numerical results in section |7.1| . 

This propagator (^) must be applied to the ground state ipo (||) in order to obtain the 
wavefunction ip{tf) at a certain time tf 

ij{tf) = U{tf,0)iJo (50) 

Applying the propagator U (t) results in just a phase factor exp{—(iE'^/h)t), and so 
the above expression (|^) simplifies to (for one period tf = 27r) 

ipix, tf) = exp( —tf) tpoix) + exp(z70) exp( —to) J U^{x, tf, y, to) tpo{y) dy (51) 

— oo 

Using the property (^ that ipoiu) approaches the spatial (5-function in the semiclassical 
limit h —>■ 0, the integral can be evaluated to 

ipix^tf) = exp( —tf)ilJo{x)+exp{i-f(f))exp{ —to)2J-U^{x,tf;0,to){52) 



19 



In order to calculate the probability amplitude p for the electron to remain bound, one 
has to project onto the ground state iPq{x) 

+00 

p= J ip{x,tf) iPq{x) dx (53) 
—00 

Using the spatial localization property (^ again, this simplifies to 

iE"^ iE"^ h 

p = exp( —tf) + exp( — to)4-f/^(0, tf, 0, to) (54) 

n n '-f 

E^ and to are known from the equations (^3]) and (^31), respectively. The phase 
is identical to because the relevant classical path x^ = ^ — cos(t) never crosses the 
origin. The first part of this expression for p clearly accounts for the background, while 
the second part determines the amplitude, phase, and period of the superposed (slow) 
modulation. 



6.2 Fundamental periodicity 

The phase of the first term in eq. ( pH ) is mainly given by the expression —Eotf/h, and 
the phase of the second term is dominated by the action S^^ /h along the classical path 
Xd = l- cos(t). 

*/ 

S'^ = J (^xli{t) + Xdit) cos(t)^dt = -z tf (55) 


The last identity is straight forward (also in atomic units). Combining the phases 
and comparing to multiples of 2 7r results (for tj = 2 7r) exactly in equation (^ with 
Ekin = 0. This is just the threshold condition for channel closing, and implies the same 
periodicity. The period in z is Az = 1/(1 + 27^), if 7 is kept fixed, or a period of 
Az = 1, if the depth of the binding potential a is kept fixed (see figure ^ later). 



6.3 More periods 



Again this result can be generalized easily to more bursts and longer final times tf > 2tt. 
Be tf = 2mT, then one has to sum over 2n bursts and amplitudes, and (this time 
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written explicitly) the result for p is 

? - 2n-l _A7 

fe=o g J2iiTh[tf -to - kn) 



— ^ 

' (to + knY + {tf — tQ — kn) cos(to) sin(to) (57) 



Ah{tf - to - kn) 



-2{to + k7T)tf + 4cos(to)(-l) - 2 + - 2 cos^(to) ) - ;^^"'(^o + kn) 



I 



Note the useful relations cos(to) = VI + 7^ and sin(to) = ^7. Propagating the system 
for longer times > 2 7r means to have more phase built up in the exponents of (pBf), 
and results in a finer resolution in the ionisation rate F; higher frequencies than the 
basic modulation period can be accounted for. How this can explain the fine-structure 
is demonstrated later on in figure 0. 
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In the above sum (^), only the contributions with k even are important. This is 
because the centers of wave packets stemming from k odd are located at about —2 to 
the left at tj = 2n7r, and therefore the overlap is very small. 



6.4 Ionization rate 

The probability w for not being ionized is now calculated by taking the absolute square 
w = and the corresponding ionization rate F, fitting the exponential decay w = 
exp{—Ttf/{27i)), can be defined as 

2n 



tf 



H\pn (58) 



For this quantity F, we will now compare the above semiclassical results with those 
from numerical simulations 
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7 Comparisons 

7.1 Comparison with numerical results 




4 6 8 10 

z at -y — O.y 

Figure 4: Numerical ionization rate (thin and jagged) and semiclassical approximation 
(thick and smooth) versus number of ponderomotive photons z for 7 = 0.7 

The numerical results are obtained by an integral equation method implemented by K. 
Sonnenmoser (|jl| for details), which allows high-resolution scans and exhibits a lot of 
fine structure. While 7 is kept fixed, z is varied, and so the semiclassical limit h 
corresponds just to 2 — 00. This means that the agreement will become the better 
the larger z is. 

The interesting region for 7 is of course 7 1, because this is the transition region 
between adiabaticity and multiphoton regime. For 7 <^ 1, the adiabaticity criterion 
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Figure 5: Numerical ionization rate (thin and jagged) and semiclassical approximation 
(thick and smooth) versus number of pondero motive photons z for 7 = 1.1 

is fulfilled, the averaged WKB- value can be justified, and is in good agreement due 
to the very construction of our theory. For 7^1, one should turn over to a pure 
multiphoton description ||23|, ^ . 



Figure ^ shows the numerical result (thin and jagged) for 7 = 0.7, together with the 
results of our theory. The background as well as the amplitude, phase, and periodicity 
of the modulation are very well comprised in the semiclassical theory for z not too 
small. 

The same is done in figure ^ for 7 = 1.1 . Here again, one recognizes that the charac- 
teristic elements of the ionization curves are in good agreement. The same is valid 

of course for all other values 7 1 and can be extended up to 7 ~ 2.5, clearly beyond 
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4 6 8 10 

z at -7 = 0.7 

Figure 6: More fine-structure by considering longer periods than 2 7r, here two cycles 
4 71 considered 

the adiabaticity regime. 

If we want to incorporate more fine-structure superposed onto the modulation, we 
can consider longer periods than 2tt. The result is a behaviour as in figure which 
resembles closely the fine-structure, even if there can be no one-to-one correspondence 
between every small wiggle. Here the wave packets stemming from ionization bursts vr 
and Svr were also taken into account, but they have very little influence on the result 
for tf = Att. The main contributions come of course from the wave packets stemming 
from t = and t = 27c. This can be seen as an example for the above quantitative 
statement about the relative importance of ionization at even or odd multiples of vr. 

The result (|56D for fixed 7 and z can easily be transcribed to other parameter combina- 
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Figure 7: Now the depth of the binding potential is kept fixed, and z is varied. The 
semiclassical theory (thick smooth curve) is compared to numerical results (thin and 
jagged). 

tions. A common representation of ionization rates F is to keep rijo = 27^2; = o? j i^Lo) 
fixed, i.e. the depth of the binding potential (last expression in atomic units as in sec- 
tion p^.l| ). One varies the intensity ^ at fixed frequency cj, which corresponds better 
to experimental situations. If one plots F over the intensity or over z = /i^/(4c(j^), one 



again obtains periodicity in the ponderomotive channel closing [^. This is depicted 
in figure |^ in comparison with numerical data; note hereby that the semiclassical limit 
no more corresponds simply to 2; — > 00, but is more involved and cannot be included 
simply into the representation. One again notices that the characteristic elements of 
ionization curves can be calculated in the semiclassical theory. 
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7.2 Comparison to other theories 

One major advantage and distinction of our theory from others (even from those claim- 
ing to be semi- or quasiclassic in some sense) is that we encounter no divergencies at 
channel closing thresholds, a result strongly supported by numerical evidence. Such 
divergencies typically occur in theories that (in our mind somehow arteficially) sepa- 
rate the ionization process into different channels, each one related to ionization by a 
distinct number of photons. Every time a channel closes, the corresponding ionization 
rate of the next higher channel (and so the overall rate) becomes infinite. 

One earlier representative (eq. [31] by Perelomov, Popov, Terent'ev |jl3|; including qua- 



siclassical features) and one more recent representative (eq. [44] by Susskind, Cowley, 
Valeo [0]; asymptotic in the number of photons for ionization) of this kind of theory 
is shown in figure |, in comparison with our result. The spikes in this figure are related 
to the closing of certain ionization channels with, say, k photons at z = Zk (eq. (0)), 
and they can be traced back to the divergence in the next higher channel with k + 1 
photons. 

In between these thresholds, the background and the imposed modulation is approx- 
imately the same. But in contrast to these more implicit theories, where results are 
quite involved, our result (|56| , ^8|) allows the direct, separate, and explicit evaluation of 



the background rate as well as of the amplitude and phase of the modulation. 
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4 6 8 10 

z at y — O.y 

Figure 8: Comparison with other theories containing divergencies at channel closing 
thresholds: semiclassical theory (thick and smooth) versus Perelomov, Popov, Ter- 
ent'ev [13] (thin hne) and Susskind, Cowley, Valeo [14] (dashed), which nearly coincide 
in this region and especially in their behaviour at thresholds. 

8 Summary 

We succeeded in a semiclassical description of time-dependent tunneling and ionization 
in an oscillating field. The characteristic features of typical ionization curves can now 
be explained using a picture of tunneling, propagating and interfering wave packets 
(formula (0)). The main ingredients are first the separation of the ionization pro- 
cess into two disctinct steps, motivated by the asymptotic evaluation of instantaneous 
WKB-rates. And second the usage of complex time and analytical continuation of 
propagators, necessary to account for tunneling by a classical path description. The 
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slow modulation with channel closing periodicity (an idea stemming from a multipho- 
ton viewpoint) can be described correctly with respect to amplitude and to phase. Even 
the rich superposed fine-structure can be accounted for by considering a multitude of 
interfering wave packets. 
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A Saddle— Point Integration Techniques 

Under the notation of saddle-point integration, there exist two different versions in the 
hterature (all considered in the limit h ^ 0): 

Method of steepest descent: 



J g{t) exp{-^f{t))dt 



^""^ -gito) exp(-l/(to)) (59) 



pito) = , ^/(to) > (60) 



Here the argument in the exponent was expanded around the minimum to of f{t), up 
to second order, and the so created gaussian integral was evaluated exactly. If there 
exist several minima, one has to take the sum over them. 

Method of stationary phase: 



J g{t) exp{^f{t))dt 



\ ^-TTT^dito) exp(^/(to)) (61) 

This is in some respect the analytic continuation of (^Sj); here every stationary point 
to with dt /(to) = is relevant, not only the minima. 

Analogy in functional integration: 

*/ 

Jvx{t)g{x{t))exp{^S[xm (62) 

u 

In this path integral, one has to integrate over all possible paths with the boundary 
conditions 

x{U) = y , x{tf) = Xf (63) 

In the semiclassical approximation, one looks for the paths Xd which make the func- 
tional S[x] stationary: 

6S[x{t)] = (64) 

One again expands around these classical paths Xd (cf. e.g. |32[) and obtains the 
so-called Van-Vleck propagator [Q, which in general must be corrected by so-called 
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Maslov factors v |^. These take into account the crossing of caustics, i.e. divergencies 

of the mixed second derivative of S'^^ . 

tt 

Q2 



= , S^' = 



V2nih\ dxfdy 
XciiU) = y , Xditf) = Xf , 6S[x(t)] 



(65) 
(66) 



B Semiclassical Propagator using the WKB— Ansatz 

This is an alternative to the construction of the semiclassical propagator f/^'^(x, t; y, ti) 
by the path intgral approach (cf. section ^]T] and appendix^. We construct U'^'^{x, t; y, ti 
in such a way that it fulfills the time-dependent Schrodinger equation (^) with respect 
to t and X up to 0{h?). 

ih^U"' = HW = ( --h^S^-h-f5(x)-xcos(t)]W (67) 
ot V 2 ox'^ } 

\YmV'\x,t-y,U) = 6{x-y) (68) 
We make the ansatz 

U = (p{x, t; y, U) exp {j^S{x, t; y, U)^ (69) 

and yield 

1 ih h"^ 

ihdt^p — ^pdtS = -{dxSy — ihOx^d^S — —cpdlS — -z-dl'^ — h'jS{x)(f — xcos{t)(f (70) 

Zi Zi Zi 

Comparing the distinct orders of h: 

h° : = dtS + ^{dxSf - X cos{t) (71) 
h' : t^S{x)if = dtV> + dx^dxS + ^^dlS (72) 
2z75(a;)^2 ^ d_^2 ^ ^{^'d^S) (73) 

The right-hand side of equation (|73|) is just a conservation equation with density ip"^ 
and flow ip'^dxS, which is fulfilled everywhere except the origin. The equation (^) for 
/i° is easily solvable, it is just the Hamilton-Jacobi equation —dtS = H{x, dxS) for the 
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motion of an electron in an oscillating electric field (with Lagrangian Lq). One can 
easily solve the appropriate equation of motion and yields the classical path Xci{t) with 
boundary conditions Xci{t) = x and Xditi) = y. 

/ \ /X / \ X — y + cos(t) — cos(ti) , . 
x,i{t) = - cos(r) + cos(t,) +y + ^ ^(r - U) (74) 

The action S is then the time integral over the classical Lagrangian Lq: 

t 

S = J Lo(x,i{t'),xa{t')^dt' (75) 



Then, as usual, the expression (fo = —d^dyS = — ti solves the equation (|72D for 
outside the origin. Now we make the ansatz 

(p = (3{x,t;y,ti)(po= (76) 

Inserting this in (^), we yield 

i-i5{x)P = dtP + 3,13 d,S (77) 

This linear partial differential equation of first order can be solved using the method 
of characteristics. The (ordinary) differential equation for the characteristic xt is 

^xt = d,S{x,t;y,ti} (78) 

But S is, as we know, the action for the classical path and therefore dxS is just the 
momentum of the classical path d^S = Xdit) (e.g. [0). So we conclude that the 



characteristic is just the classical path xt = Xd- We obtain the following (ordinary) 
differential equation for (3 : 



d 



dt' 



d 



-J = dt(3 + d,f3 -Xd = i^Six) (3 (79) 



dt 



This can be integrated straight forward and the result is 

t 

(3{t)l = exp («7 / 6{xd{t')) dt') (3^ (80) 



t 
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This again can be understood as a phase jump for every crossing of the 5-potential 



at the origin; the constant Po is fixed to /3o = 1/ V^mh by the required normahzation 
of U'^. 



To summarize, we end up with the following result, identical to eq. (^) 

1 i * 

U'%x, t; y, U) = , exp (^S[xci] (x, t; y, U)) exp (i-f f 6{xci{t')) dt') (81) 

C Complex Time Description of Tunneling Processes 



The usefulness of a complex time coordinate for describing tunneling processes is best 
demonstrated by a simple example. Consider an inverse harmonic potential barrier 
V{x) = — x^/2 and a particle coming in from — oo with energy E = —1/2. The classical 
equation of motion Xci= Xd is fulfilled by the classical trajectory Xci{t) = — cosh(t), 
restricting the particle to Xd < —I- 

In this time-independent problem, the total energy £" is a constant of motion: 

E=^--x^ = -- (82) 
2 2 2 ^ ^ 

This differential equation can be solved for t{x): 

f dx' 

t{x) = ± / ^^^== (83) 



-1 



V2^ 



X' 



The range of x can now be extended formally to x > — 1, then t acquires an imaginary 
part for x' G [—1,1] 



dx' 



-1 



V2ET 



X' 



/2 



in (84) 



The tunneling trajectory through the barrier can be described in the complex t-plane 
by the path [—oo, 0], [0, in], [i n, +oo + in]. The second part is the tunneling process 
through the barrier, noting that 

XdiiT) = — cosh(ir) = — cos(r) 

The third part describes the propagation on the other side of the barrier, noting that 

Xdit + in) = — cosh(t + in) = + cosh(t) 
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motion in x-space 



complex t-plane 




X 



Im(t) 



iPi ^ 



Tunneling 







Re(t) 



Figure 9: Tunneling in x-space and the appropriate complex time path for an inverse 
harmonic potential 

Figure ^ illustrates this by contrasting the path in x-space to the path in the complex 
t-plane. For mathematical foundations of the analytic continuation applied above see 



for recent applications of this method compare e.g. [36, 42 
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D Extensibility and Comparison to other Models 

In section there was explicitly made use of the scaling property 6{Px) = 5{x)/[3 
(for /5 > 0) of the one-dimensional 5-function. For the Coulomb-potential Vc, the 
scaling behaviour is identical, so that the whole calculation can be repeated. The 
main difference consists in a different phase factor ipc (cf. (|32|)), which is now the line 
integral over the Coulomb-potential: 

*/ 

(t>c = j VcMt))dt (85) 

u 

This is due to the fact that the binding potential is suppressed by a factor h in the 
transformed Schrodinger equation (^). As derived in sections and appendix |B|, such 
a suitably suppressed binding potential does not influence the classical trajectories, 
but only changes the phase transported along these trajectories by an additional phase 
factor (pc (given above). 

This separation is no longer valid for other types of potentials with different (or without) 
scaling properties. Here the classical equations of motion must be solved fully for 
the binding potential plus electromagnetic field. In the time-dependent case, there 
generally does not exist a first integral of motion (like the energy in the static case), 
and so the classical system is not integrable in closed form. 

The physics of the ionization process should remain the same: complex time-dependent 
tunneling followed by free propagation (two-step models). This is reflected by the fact 
that ionization curves for other model potentials show qualitatively the same charac- 
teristics. The semiclassical "sum over classical paths" -method claims to work in this 
case anyway, but unfortunately this cannot be done in analytical form because of the 
nonintegrability already on the classical level. 
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